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ABSTRACT 

A formula for the mass-gap of the supersymmetric CP n_1 sigma model (n > 1) in 
two dimensions is derived: m/A^r = sin(7rA)/(7rA) where A = 1/n and m is the mass of 
the fundamental particle multiplet. This result is obtained by comparing two expressions 
for the free-energy density in the presence of a coupling to a conserved charge; one expres- 
sion is computed from the exact S-matrix of Koberle and Kurak via the thermodynamic 
Bethe ansatz and the other is computed using conventional perturbation theory. These 
calculations provide a stringent test of the S-matrix, showing that it correctly reproduces 
the universal part of the beta-function and resolving the problem of CDD ambiguities. 
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1. Introduction 



This paper is the second of a pair (see [1]) concerned with exact results for super- 
symmetric sigma-models in two dimensions. The general strategy, which has already been 
successfully applied to non-supersymmetric theories in [2,3,4,5,6,7,8], is to calculate the 
free energy of a theory in the presence of a background field by using the Thermodynamic 
Bethe Ansatz (TBA) equations derived from some postulated exact S-matrix. By com- 
paring this result with a calculation of the same quantity in standard perturbation theory, 
the validity of the S-matrix can be checked and the exact mass-gap for the model can be 
extracted. In [1] the supersymmetric 0(A) sigma model (N > 4) was considered and it 
was shown how the novel difficulties associated with a supersymmetric theory — in partic- 
ular the problem of diagonalizing the resulting TBA system — could be overcome for this 
simplest family of examples. In this sequel we consider the supersymmetric CP n_1 sigma- 
models. These theories differ from the 0(A) models in having A = 2 supersymmetry 
and also in a number of other important respects. For our purposes, the most important 
difference is that the TBA equations we shall have to solve involve anti-particles as well 
as particles. We shall see that the methods described above can, nevertheless, be applied 
successfully in these cases too. 

The bosonic CP n_1 models (n > 1) first attracted attention as "toy" field theories 
with instanton solutions [9] generalizing the bosonic 0(3) model (the 0(3) and CP 1 models 
coincide but there are no instanton solutions in the 0(A) models for A > 3). The analogy 
with QCD is reinforced by the fact that the bosonic CP n_1 theories are asymptotically 
free, they generate their mass dynamically, and they are confining. Unfortunately, the 
classical integrability of the bosonic CP n_1 models appears to be vitiated at the quantum 
level by anomalies in the conserved currents. For the supersymmetric CP n_1 models 
[10,11,12] on the other hand, integrability is maintained at the quantum level due to a 
cancelling of the anomalies between the bosonic and fermionic degrees of freedom [13] 
and exact S- matrices have been proposed by Koberle and Kurak [14] . The structure of the 
quantum supersymmetric CP n_1 models is simpler than that of their bosonic counterparts 
in that they do not display confinement [10,11]. Even so they are highly non-trivial, 
asymptotically-free quantum field theories with dynamically generated mass, and they 
have been used — along with their bosonic counterparts — to gain some profound insights 
into non-perturbative effects in field theories, such as the relationship between instantons 
and the 1/n expansion (see eg. [9,10,11,15,16] and references therein.) More recently they 
have also attracted attention as massive integrable perturbations of A = 2 superconformal 
field theories and as useful examples to which new non-perturbative techniques can be 
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applied [17,18,19]. 



The model and its S-matrix 



The supersymmetric CP n 1 model is defined by a lagrangian density 
1 

2^ 



1 

+ 4 



where A M = i «<9^ a - ^a^a) • 



(2.1) 



We work throughout in Minkowski space and our conventions agree with those of [20] . The 
fields z a and ip a are an n-component complex scalar field and an n-component complex 
Dirac fermion, respectively, which satisfy the constraints z*z a = 1 and z*tf) a = 0. The 
theory is clearly invariant under global SU(n) transformations on these fields. There is 
also a local U(l) invariance under which the quantity above transforms as a gauge field; 
this means that the theory can be interpreted as a sigma-model whose target manifold is 
the complex projective space CP n_1 . 

Applying well-known general results (see eg. [20] and [21]), the two-loop beta-function 
for this model and the corresponding behaviour of the running coupling constant can be 
written 

[3(g) = -Pi9 2 - fog 3 + 0(g 4 ), 

g(n/A) A ft A V ln(/VA) / 

where (3i = n/n, ft = 0. 

From these expressions we see that the theory is asymptotically free, with dynamical mass 
generation. 

It is natural to expect the spectrum of the theory to contain supersymmetric multiplets 
of particles, some of which are "fundamental" , in the sense that they carry the quantum 
numbers of the fields in the classical lagrangian, and others which can be regarded as 
bound states. We denote the quantum states corresponding to the fundamental particles 
by |a,z,6>), where i = 0, 1 denotes a boson or fermion respectively, a is the SU(n) vector 
index of the n dimensional representation and 9 is the rapidity of the particle, so its 
velocity is v = tanh(6'). The spectrum also includes a set of fundamental anti-particles 
\a, i, 6) transforming in the n representation of SU(n). It turns out that these fundamental 
anti-particles can be formed as bound-states of the fundamental particles (or vice-versa; 
an example of "nuclear democracy"). 
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The integrability of the model [13] implies that the S-matrix factorizes and that all S- 
matrix elements can be deduced from the two-body one [22]. On the basis of this, Koberle 
and Kurak [14] (see also [17,23,24,25]) proposed an S-matrix to describe the scattering of 
the fundamental particle multiplet in accordance with the SU(n) symmetry of the model 
and all the usual axioms of S'-matrix theory. Their proposal can be written in an illumi- 
nating way in which the supersymmetric and SU(n) degrees of freedom are factored. This 
was made explicit in [17] following the discussion for a general supersymmetric theory in 
[26]. 

In detail, the two-body S'-matrix elements for the fundamental particles can be written 



(c, k, 9 2] d, /, ft, out|a, i, i; 6, j, 2 , in) = SW =2 (0i - 2 ) • 7 'S C gn(0i - 6 2 ) c a d b 



\cd 



(2.3) 



The SU(n) part of the S-matrix is the factorizable S-matrix of the fundamental vector 
particle of the SU(n) chiral Gross- Neveu model [27]: 



Scgn(0)S6 = ^i(0) 



5 ad 5 bc 



2tuA 

e 



5 ac 5 bd 



(2.4) 



with A = 1/n and with the unitarizing/crossing scalar factor 

r(l + i0/27r)r(-A - i0/27r) 



Yi(9) = 



r(-z0/2yr)r(l- A + i0/27r)' 



(2.5) 



The supersymmetric part of the S-matrix is the N = 2 supersymmetric Z n minimal S- 
matrix [17] which has the form 



S N =2(6) = Y 2 (6) 

/sinh(0/2 + Z7rA) 



V 



x 





zsin(7rA) sinh(0/2) 

sinh(0/2) zsin(TrA) 











\ 



(2.6) 



sinh(0/2 - in A) / 

where the rows and columns are labelled in the order (0, 0), (0, 1), (1, 0), (1, 1) and where 

1 



Y *W) = - 



sinh(0/2 + Z7rA) 
T 2 (i6/2n + j)T(-i6/2n +j + A)T(-i6/2n + j - A) 

X 11 



(2.7) 



=i r 2 (-z0/2yr + j)T(i6 /2tt + j + A)T(ie/27r + j - A) 



The S-matrix elements of the fundamental anti-particles can be found by crossing the 
elements (2.3). 

It is important that in (2.3) we have chosen an ordering in the final state where the 
particle of rapidity 9 2 is to the left of the particle of rapidity 6>i; it is only this "modified" 
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S- matrix (using the nomenclature of [26]) which exhibits the factorization between the 
supersymmetric and bosonic degrees of freedom manifested in (2.3). This is opposite to 
the choice made in the original treatment of [14], which means that some care is required 
in comparing the signs for amplitudes involving two fermions. It is also worth pointing out 
that although it is this same "modified" S-matrix which appears in [17], we are including 
the simple pole in the SU(n) factor rather than in the supersymmetric factor; the net result 
is easily seen to agree with [17]. 

The expression (2.3) is "minimal" in the sense that it has the minimum number of 
poles and zeros on the physical strip (the region < Im(9) < n) consistent with the 
requirements of symmetry, existence of a bound-state and the axioms of S-matrix theory. 
But this still leaves open the possibility of adding CDD factors to the S-matrix; these 
spoil none of the axioms, they introduce no new poles on the physical strip and they 
passively respect the bootstrap equations. For our model the CDD ambiguities correspond 
to multiplying the S-matrix of the fundamental particles (2.3) by factors of the form 



where < a < 2. One of the conclusions of this paper will be that the minimal form (2.3) is 
the true S-matrix of the theory, so that all CDD factors are ruled out. As in previous work 
[1-8] this will follow from the consistency of our calculation using the minimal S-matrix 
with a calculation in perturbation theory. 

Although it will not concern us directly here, we point out that the complete spec- 
trum of the model can be determined using the bootstrap procedure. The S-matrix (2.3) 
for the fundamental particles has a simple pole on the physical strip at 9 = 2-7rzA which 
corresponds to a bound state transforming in the antisymmetric tensor representation of 
SU(n). Continuing the bootstrap in this way one finds a spectrum of bound-states which 
is identical to the SU(n) chiral Gross-Neveu model, namely m r = msin(7rrA)/ sin(7rA), 
1 < r < n, where the r th bound-state transforms in the r th fundamental representation 
of SU(n) (each particle carries in addition supersymmetric quantum numbers). The fun- 
damental anti-particles correspond to the n— 1 th bound-state. We shall only require the 
S-matrix elements of the fundamental particles and the fundamental anti-particles for our 
calculation. 

3. Coupling to a conserved charge 

As summarized in [1], we wish to couple the model to a chemical potential h via a 
conserved charge Q and to calculate the resulting free energy as a function of h. In other 




(2.8) 



sinh(f + g(2-a))sinh(§ 
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words, we seek the response Sf(h) = f(h) — /(0) of the ground state energy density for 
the system with Hamiltonian modified from H to H — hQ. The TBA equations allow 
one to find an expansion for 8f(h) of the form h 2 Fi(h/m) valid when li> m and in this 
regime one can calculate the same quantity in perturbation theory to obtain an expression 
/i 2 i ? 2(/i/A), where A is the usual dimension-full parameter of perturbation theory. Setting 
F\(h/m) = F2(h/A) and comparing the first few leading order terms gives a powerful test 
of the S-matrix and allows one to extract the mass gap m/A. 

The analysis of the TBA system is only tractable if we can choose Q so that the 
new ground-state consists of a restricted number of particle types. In [1] we found that 
the presence of supersymmetry complicated matters because the ground state formed a 
supersymmetric doublet with non-trivial scattering, unlike the examples in [3-5] where one 
could select Q so as to give a single particle type. Despite this complication, the problem 
considered in [1] could still be solved. In the present case we might be tempted to make 
the same choice as for the SU(n) principle chiral model [4] and the chiral Gross-Neveu 
model [6]: 

Q = diag(l, , (3.1) 

\ Ti — 1 n — 1 J 

for which the fundamental multiplet \l,j,9) has the largest charge/mass ratio. As we 
shall see in the next section, however, this choice is inconvenient because, even if the TBA 
system proved tractable, it would require a three-loop perturbation theory computation 
to extract the mass-gap. To avoid this we are motivated to consider the alternative choice 

Q = diag(l,-1,0,...,0). (3.2) 

With this coupling the situation is still more complicated than in [1] because there are now 
two fundamental doublets with the largest charge/mass ratio, namely \l,j,9) and \2,j, 9), 
a feature which clearly arises because of the presence of distinct antiparticles. Nevertheless 
we shall see that the resulting TBA system can be analyzed successfully and that the mass 
gap can again be extracted by comparison with a perturbation theory calculation to just 
one loop. 

It is worth re-iterating the point already emphasized in [1] that it is an assumption 
that only 9) and |2, j, 9) appear in the ground state and that in particular no bound 
states contribute. As in most previous work of this type, we rely on the consistency of 
our final results to vindicate this assumption. An important property of the states which 
we are assuming appear, is that their scattering is purely elastic in the space of SU(n) 
quantum numbers. Of course the scattering is still non-diagonal in the supersymmetric 
subspace. The SU(n) part of the S-matrix elements of 1 with 1 and 2 with 2 coincide: 

q m ii q (f)] -22 r(l + zfl/27r)r(l-A-zfl/27r) 

ScgMu = SccAOh-2 = vK^W^W^VWm- (3 - 3) 
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The SU(n) part of the S-matrix for 1 interacting with 2 can be found by crossing: 

<? r^i_c r(9 u2_ rq-«g/27r)ru-A + «g/27r) 

5 C GN((?) ia - ^CGN(^) 51 - r(i+ ^ /27r)r(i _ A _^ /27r) - ( 3 - 4 ) 

A crucial point is that there is no reflection amplitude for the interaction of 1 and 2: the 
scattering is purely elastic. 



4. Free-energy from perturbation theory 



To couple the theory to the charge (3.2) by changing the Hamiltonian H — > H — hQ 
we can make the replacement do — > do + ihQ in the Lagrangian (2.1). It will turn out to be 
sufficient to calculate the ground-state energy density of this theory to one loop. We must 
therefore expand the Lagrangian to quadratic order in an independent set of fields and we 
can drop all terms which are independent of h to this order because we are interested only 
in the response of the free-energy density, 5f(h) = f(h)—f(0). By exploiting the local U(l) 
invariance of the action, we can take z\ to be real and we can solve the bosonic constraint 
by writing z\ = a/(1 - |-7T | 2 ) ( i + <j>) and z 2 = e 10 ' yj '(1 - |7r| 2 )(i - 0) where 7r = (z 3 , ...,z n ) 
and 9, (j) are real. The fermionic degrees of freedom and the variable 9 decouple to this 
order, and we are left with the expression 



'1— loop 



^ {{d^f + |<V| 2 + h 2 - 4/iV - h 2 \n\ 2 } 



(4.1) 



Using dimensional regularization with the MS-scheme gives the one-loop free energy 

■(1-Hh 2 /^ 2 )). (4.2) 



h 2 h 2 , „ nh 2 



^ = -2^-V ln2+ 4. 

We now substitute for the running coupling. The result, expressed in terms of the one-loop 
and two-loop beta functions coefficients, is: 



Sf(h) 



■A 



In 



h 



12 In 2 ^ h 

~ n H 15~ + ~m ln ln ~\ 1 - ° 

2 7T/?i (3f A m 



ln ln(h/A 



\n(h/A- 



MS/ 



(4.3) 



Substituting the specific values for these coefficients given in (2.2) we obtain 



Sf(h) = 



h 2 n 



ln 



h 



1 21n2 _ 

- + + o 

2 n 



Hh/Aus) 



(4.4) 



Notice that if we chose the charge Q to be (3.1) then there would be no "tree-level" 0(1/ g) 
term. In this case it would be necessary to do a three-loop computation in order to extract 
the mass-gap. 
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5. 



Free-energy from the S-matrix 



In this section we write down the TBA equations for the model and solve them in the 
limit h~^> m. We must follow the hypothesis made earlier that only the multiplets |1, j, 9) 
and 1 2, j, 9) contribute to the ground-state. Since the scattering of these multiplets is purely 
elastic, it will not be necessary to perform a diagonalization in the space of SU(n) quantum 
numbers (although this diagonalization can be done [28]). The remaining difficulty is that 
the S-matrix for these favoured states is still non-diagonal in the supersymmetric subspace. 
But this problem can also be solved since it has been shown by Fendley and Intriligator 
[17] that it is equivalent to diagonalizing the transfer matrix of the six vertex model at the 
free fermion point. 

In [17] a system of equations is derived which links the density of single particle 
states in rapidity space g a (9) to the density of occupied single particle states o" a (6>), where 
a = 1, 2. The equations also involve densities Pj*~(9) and P l ~(9) : I = 0, corresponding to 
two "supersymmetric magnons" , reflecting the fact that there are two supersymmetries: 

Qa(0) = —cosh 6 + <j) ab *a b (6) + (f> a i * P+ (0), / x 

2tt (5.1) 

p+{9) + p l -{e) = <t> a i*a a {e), 

where / * g{9) = f_ d9' f{9 — 9')g{9'). The kernels appearing here are 

^) = ^^lnScGN(«, 

<M0) = fco(0) = ° in(7rA) , AV (5.2) 
coshfc' — cos(7rA) v ' 

a <a\ a (o\ sin(TrA) 



cosh 9 + cos(7rA) 



where as before a, b = 1,2 and the elements Scgn^)^ are written down in (3.3) and 
(3.4). Notice that the kernel which multiplies the densities of the particles a a {9) involves a 
contribution only from the chiral Gross-Neveu part of the S-matrix — so if we removed the 
magnon terms the equations would be identical to those for the chiral Gross-Neveu model. 
Notice also that the super-magnons do not interact amongst themselves. 

To find the TBA equations [29] we define the "excitation energies" of the particles 
e a {9) and the magnons £i(9) at finite temperature T by 

Qa{9) e«-W/ r + l' P+(9) ' 1 } 
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We shall require the TBA equations at zero temperature with the field h coupling via Q 
acting as a chemical potential. At T = 0, e a (±6>£) = 0, where 6>p are Fermi rapidities and 
e a (9) is negative precisely when — 0% < < 6>p. It is convenient to introduce the following 
notation 

/±(6>) = (/(*) m<° (5.4) 

1 otherwise. 

The free-energy per unit volume at T = is 

/oo 
d0[e^(0) + e^(0)]cosh0, (5.5) 



where e a (6>), a = 1, 2, are the solutions of the T = TBA equations: 

fa (0) - <Pab * ^ (°) ~ <f>al * f j" (^) = 771 COsh 61 ~ ^' 

&(0)-0a,*eJ(0) = O. 



(5.6) 



To simplify (5.6) it is important to notice that (f> a i(0) is a positive kernel, which implies 
that the magnon variables are given by £ z + (0) = and £f (t) = (p a i*^a (#) • Furthermore, the 
solution does not distinguish between the values of the favoured SU(n) quantum numbers 
and so we have ci(6) — e§(#) = e(0). The four equations in (5.6) then reduce to a single 
equation for e(6>): 

e + (0)+R*e-(0) = mcosh0-h, (5.7) 

where the kernel is 

R(0) = 6(0) - M0) - l2 (0) - [0io + 0io] * [0io + 0io] (5-8) 
and the expression (5.5) for the free-energy density becomes 

/oo 
d0e _ (0)cosh0. (5.9) 
-oo 

The Fourier transform of the kernel in (5.7) is 

m = f°° *± co S (^) C ° sh " 1 - 2A 'r /2)Sinh( " A " > e^, (5.10) 
Jo ^ cosh (ttoj/2) 

where A = 1/n. As in [1] we see that this Fourier transform vanishes at the origin so 
that the solution resembles those for the bosonic models discussed in [3-4]. We now seek 
an expression for the Fourier transform of the kernel 1/(G+(uj)G-(uj)) where G±(ui) are 
analytic in the upper (lower) half planes and G-(u) = G + (—uj). The unique solution is 

r ( , A ~ »(1 ~ 2A)o;/2)r(l - iAu) i ln( _ lAoj) 

+ {) (5.11) 

x -iw(^-A)(l-ln(-iw(^-A)))-iwA(l-ln(-iwA))+ia;(l-ln(-i|w)) 
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From [4] we know that if G + (i£) has an expansion for small £ like 



G + (i0 = ^e- a ^(l-bZ + O(e)), 
then the free-energy density for h ^> m takes the form 



Sf(h) = 



2u2 



h 2 k 



, h i 

In h In 

m 



'2nke~ 



G+(i) 



h 



l + a(7E-l + ln8) 
lnln(/i/m) 



+ (a+ i)lnln— + . , 

m \ m(/i/m 



Our kernel does indeed have an expansion of the form (5.12) with 
k 



1 1 V2nke~ b sin(TrA) 7e/2+ ( |+2a) in 2 

G+(i) 



/ttA' 2 
The resulting expression for the free energy is 



TTA 



J v ' 2ttA [ m V ^ A / 2 V hi(/i/m) 



(5.12) 



(5.13) 



(5.14) 



(5.15) 



6. 



Comparison and conclusions 



Comparing (5.15) with (4.4) we see that the result from the TBA calculation correctly 
reproduces the universal coefficients of the beta-function and we extract the following value 
for the mass-gap for the supersymmetric CP n_1 model 



m 



sin(7rA) 

7rA 



n 



n > 1. 



(6.1) 



It is instructive to consider what would have happened if we had chosen the other 
charge (3.1), rather than (3.2). In that case only the multiplet 9) would have appeared 
in the ground-state and the resulting TBA equation would have had the same form as (5.7), 
but with a different kernel: 



R{9) = 5(9) - M9) - 1O * 0io(0) - 0io * <Pro(0). 



(6.2) 



The Fourier transform of this kernel does not vanish at the origin, and so the situation 
is analogous to the Gross- Neveu models [6,5]. In these cases, the expansion of the free- 
energy density for h ^> m has a different form which has a well-defined limit as h — > 00, 
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unlike (5.15). This matches precisely the fact that with this different choice of charge the 
perturbative expansion of the free-energy density is also markedly different since there is 
no "tree-level" contribution and it would require a three-loop calculation to extract the 
mass-gap. 

In addition to our result for the mass gap, an important conclusion of our paper is that 
we have resolved the problem of CDD ambiguities in the ansatz of Koberle and Kurak. 
Any additional CDD factors of the form (2.8) would alter the kernel R(0) of the TBA 
equation and the thermodynamics of the system depends so sensitively on this that we can 
argue with some confidence that the solution for the free-energy density would not match 
that found in perturbation theory. 

It is interesting to compare our result for the mass-gap with the spectrum for the 
super-CP n_1 model proposed in [18] on the basis of very different methods. It is clear 
that the the n-dependence of these expressions is in complete agreement (setting r = 1 
in the formula given in the discussion following equation (6) in the first reference in [18] 
should give the mass of the fundamental multiplet) but the connection between the overall 
scale factors is less clear. It would be interesting to investigate in detail how the results 
of [18] could be expressed in terms of the more conventional scale parameter used in this 
paper. 

J.M.E. is grateful to Michele Bourdeau for discussions concerning [17-19]. 
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